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ABSTRACT 

Although  there  has  been  interest  in  the  behavior  of  metal  plates  under  blast  and  projectile 
loading  for  many  years,  definitive  open-source  analysis  has  only  been  recently  forthcoming. 
This  analysis  is  most  often  in  the  form  of  scaled  recreations  of  the  dynamic  blast  event,  or 
“live  fire”  tests.  New  developments  in  methods  of  recreating  blast  and  projectile  induced 
plate  failure  using  a  quasi-static  approach  provide  possible,  accurate,  alternatives  to  the 
cumbersome  and  expensive  live  fire  test. 

This  research  endeavors  to  develop  an  accurate,  quasi-static  method  of  recreating  the 
petalling  phase  of  blast  and  projectile  failure  in  metal  sheets,  based  on  a  modified  trousers- 
type  test.  By  using  the  trousers-type  fracture  test  the  overall  plastic  bending  kinematics  of 
the  fractured  petal  is  preserved,  as  well  as  the  mixed  mode  (mode  one  and  mode  three) 
fracture. 

Through  analytical  and  qualitative  analysis,  a  testing  apparatus  to  generate  this  trousers-type, 
plastic  bending  and  mixed  mode  fracture  was  designed  and  machined. 
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Rh 

R-mm 

Rmax 

T 

Wb 

Wm 

wt 
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Rectangular  tab  width. 

Parallel  pre-cut  length. 

Crack  tip  opening  angle. 

Crack  tip  opening  displacement. 

Total  instantaneous  force  exerted  in  one  petal. 

Wedge  flap  bending  force. 

Wedge  cutting  friction  force. 

Wedge  cutting  membrane  force. 

Total  instantaneous  force  exerted  in  one  rectangular  tab. 

Total  minimum  wedge  cutting  force. 

Minimum  instantaneous  wedge  cutting  force. 

Panel  geometry  parameter  from  Office  of  Naval  Research  (ONR)  damage 
prediction  model. 

Plate  thickness. 

Instantaneous  length  of  petal  hinge  line. 

Material  properties  parameter  from  ONR  prediction  model. 

Fully  plastic  bending  moment  per  unit  length. 

Number  of  symmetric  petals  in  general  petalling  geometry. 

Resultant  hole  size  from  ONR  damage  prediction  model. 

Minimum  predicted  hole  size  from  ONR  damage  prediction  model. 
Maximum  predicted  hole  size  from  ONR  damage  prediction  model. 

Plate  thickness  from  ONR  damage  prediction  model. 

Bending  work  dissipated  in  one  petal. 

Membrane  work  dissipated  in  one  petal. 

Total  work  dissipated  in  one  petal. 

Total  minimum  wedge  cutting  work. 

Distance  from  instantaneous  crack  tip  along  crack/ fracture. 

Instantaneous  length  of  plastic  zone  near  crack  tip. 

Angle  of  crack/ fracture  convergence. 
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Instantaneous  local  crack  width. 

Instantaneous  crack  tip  opening  distance. 

Non-dimensional  CTOD  parameter. 

Crack  tip  opening  displacement  (CTOD)  parameter. 

Cross  head  vertical  displacement. 

Cross  head  vertical  speed. 

Plastic  bending  moment  amplification  factor. 

Central  petal  semi-angle  in  general  petalling  geometry. 
Cutting  wedge  semi-angle. 

Instantaneous  length  of  crack  or  fracture. 

Instantaneous  length  of  petal. 

Instantaneous  petal  length  rate  of  change. 

Pre-cut  petal  length. 

Instantaneous  radius  of  curvature  of  petal  at  the  hinge  line. 
Rolling'cylinder  inner  radius. 

Rolling  cylinder  outer  radius. 

Average  flow  stress. 

Instantaneous  rotation  of  petal  at  hinge  line. 
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INTRODUCTION 

It  is  inherent  in  the  design  of  any  warship  to  provide  robust  resistance  to  hull  and  ship 
system  damage  under  batde-type  conditions.  Since  the  extensive  naval  engagements  of 
World  War  II  there  has  been  a  sustained  effort  to  study  the  detailed  batde  damage  reports 
of  naval  vessels  in  the  Pacific  Theatre  with  the  goal  of  understanding  the  mechanics  of  their 
damage  and  failure.  This  analysis  led  to  the  development  of  many  protection  systems,  to 
abate  the  damage  inflicted  by  gunfire,  torpedo  and  mine  attack.  But  as  naval  weapon 
technology  rapidly  developed  in  the  post-World  War  II  years,  into  the  Cold  War  era  and 
beyond,  the  damage  mitigation  systems  have  not  kept  pace.  Little  is  known  of  the  effects 
of  modem  naval  weapons,  such  as  anti-ship  cruise  missiles,  advanced  capability  torpedoes, 
and  shaped  charge  warheads,  beyond  the  largely  classified  data  provided  by  full-scale 
weapons  tests  on  obsolete  platforms.  Even  less  is  known  about  the  batdefield  efficacy  of 
the  modem  systems  designed  to  counter  these  new  weapons. 

The  most  recent  data  point  for  analysis  is  the  damage  of  the  U.S.S.  COLE  (DDG-67) 
on  12  October  2000  in  the  port  of  Aden  Yemen.  It  is  unofficially  estimated  that  the  state- 
of-the-art  Arleigh  Burke-class  Guided  Missile  Destroyer  was  rocked  by  between  400  and 
700  pounds  of  C-4  explosive  detonated  at  the  waterline,  at  a  standoff  of  0  to  10  feet  from 
the  hull.  The  extent  of  the  damage  to  the  ship  can  be  clearly  seen  in  Figure  1,  showing  the 
20-foot  by  40-foot  hole  tom  into  the  port  side  hull  of  the  ship. 

As  can  be  seen  in  this  figure,  a  preponderance  of  the  damage  occurred  below  the 
waterline,  and  the  overall  characteristics  of  the  damaged  area  were  similar  to  the  findings  of 
Cole  (1948,  [24]),  Wierzbicki,  et  al.  (1996,  [18]  and  1999,  [3]).  The  blast  resulted  in  a 
spherical  bulging,  or  dishing,  of  the  hull  plate  prior  to  the  onset  of  tearing,  or  petalling. 
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Figure  1:  USS  Cole  Port  Side  Damage  (from  U.S. 

Navy  Information  Office) 

The  USS  COLE  was  designed  using  the  U.S.  Navy  survivability  standards  set  forth  in  a 
series  of  Design  Data  Sheets  (DDS’s),  specifically  DDS  079-1  (1976)  “Stability  and 
Buoyancy  of  U.S.  Navy  Surface  Ships,”  DDS  072-3  (1988)  “Conventional  Weapons 
Protection  (fragments),”  DDS  072-4  (1986)  “Hull,  Mechanical,  and  Electrical  Systems 
Survivability,”  DDS  072-6  (1987)  “Shaped  Charge  Warhead  Weapon  Effects  Data,”  DDS 
072-7  (1988)  “Conventional  Airblast  (proximity),”  and  DDS  072-8  (1986)  “Conventional 
Airblast  (contact  and  internal)  Design  and  Analysis  Methodology.”  These  design  guidelines 
undertake  to  outline  a  systems-based  approach  to  the  mitigation  of  damage.  They  were 
conceived  using  classified  explosive  deformation  and  holing  studies  in  naval  vessels, 
empirically  based  on  data  accrued  through  years  of  live  fire  tests  conducted  by  the  Office  of 
Naval  Research  (ONR). 


Figure  2:  Current  Stiffened  Panel  Damage  Prediction 
Model 
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The  resultant  empirical  engineering  tool  developed  by  ONR  (Figure  2)  suggests  the 
general  relationship: 


Rmin*Rh*Rmax=  HGJM) 


Where: 

Rh  =  Resultant  Hole  Size 

R^in  =  Minimum  Predicted  Hole  Size 

R^  -  Maximum  Predicted  Hole  Size 

G  =  Panel  Geometry 

T  —  Plate  Thickness 

M  =  Material  Properties 


The  direction  of  this  study  is  to  bring  further  illumination  to  the  characteristics  of  T 
and  M,  Material  Properties,  in  the  above  relationship.  This  research  is  primarily  concerned 
with  the  cracking  and  petalling  phase  of  fracture  of  hull  plating  subjected  to  a  contact, 
underwater  or  air  explosion.  It  will  serve  to  augment  previous  work  in  relating  blast-type 
failure  of  metal  plate  using  a  quasi-static  approach.  The  objective  is  to  provide  a  method  to 
more  easily  obtain  accurate  data  on  the  material  properties  of  steel  plate  for  this  mode  of 
failure. 


STATEMENT  OF  PROBLEM 

The  investigation  of  holing  failures  in  naval  plate  steel  has  been  ongoing  since  the 
transition  from  wooden  ships  to  steel,  around  the  turn  of  the  last  century.  The  basis  of 
most  research  in  the  field  began  with  the  goal  of  protecting  naval  ships  from  the 
penetration  of  artillery  shells.  Early  research,  conducted  by  Bertram  Hopkins  (1912), 
examined  the  resistance  of  various  armor  plating  to  ballistic  particle  penetration.  His 
findings  were  among  the  first  to  illustrate  the  geometry  of  holing  failure  in  metal  plates, 
including  plate  dishing,  and  petalling  from  the  formation  of  radial  cracks.  Figure  3. 


Figure  3:  Armor  Plate  with  Artillery  Penetration  (from 
Atkins  et  al.  [15]) 

With  the  experiences  of  the  two  World  Wars  came  concern  for  holing  failure  in  naval 
ships  from  the  explosive  force  of  torpedo  attack.  Taylor  (1948  [21])  and  Cole  (1948  [24]) 
conducted  a  comprehensive  study  of  submerged  blast  waves  and  their  effects  on  thin  plates 
that  formed  the  analytical  basis  of  all  current  blast  damage  prediction  methods.  Although 
most  subsequent  research  into  this  field,  conducted  by  ONR,  has  been  classified 
confidential,  open-source  study  has  been  conducted  on  plate  tearing  and  petalling  caused  by 
on-contact  explosives  by  Keil  (1956  [31]  and  1961  [32]),  Nurick  (1996  [20]),  Wierzbicki 
(1996  [18]  and  1999  [3]),  and  Rajendran  et  al.  (2001  [1]).  The  most  comprehensive  research 


program  in  the  perforation  of  plates  by  projectiles  was  conducted  by  Goldsmith  et  al.  (1978 
[38],  1983  [34],  1984  [35],  1984  [35],  1984  [36],  1984  [37]). 

Through  this  not  insubstantial  body  of  data,  the  characteristics  of  mild  and  high 
strength  steels  have  been  extensively  documented;  however  no  simple,  reliable  method  of 
predicting  hull  plate  blast  damage  has  been  developed.  Although  computer  codes  for  the 
prediction  of  blast  damage  are  available,  none  provide  more  than  a  rough  estimate  of 
potential  damage.  As  a  result,  nearly  all  of  the  definitive  blast  damage  prediction  is 
conducted  using  scaled,  live  fire  tests,  requiring  substantial  time  and  resources. 

Within  the  last  few  decades  there  has  been  a  drive  to  characterize  and  study  the  effects 
of  these  dynamic  failure  events  using  a  quasi-static  approach.  This  quasi-static  approach  to 
the  issues  of  ballistic  penetration  and  blast  failure  of  metal  plates  has  two  purposes: 

1.  To  relate  the  time-pressure  history  of  the  dynamic  event  to  the  corresponding 
force-displacement  history  of  the  quasi-static,  and  in  so  doing  relate  the 
incident  blast  wave  energy  direcdy  to  the  plastic  deformation  and  fracture  in 
the  material. 

2.  To  work  toward  development  of  a  fundamental  crack  propagation  criterion 
through  the  examination  of  crack  initiation  and  propagation  and 
corresponding  incremental  strains. 

These  two  purposes  work  toward  the  goal  of  improving  existing  computer  finite  element 
codes,  leading  to  improved,  simplified  and  reliable  damage  prediction  tools. 

To  that  end,  research  has  been  conducted  relating  the  ballistic  particle  holing  failure 
mode  using  a  quasi-static  method.  Most  recently,  Atkins  (1998  [15])  used  conical  and 
spherical  penetrators  to  observe  the  necking,  initial  fracture  (disking),  and  radial  cracking 
(petalling)  in  ductile  materials.  Amdt  et  al.  (2001  [16])  conducted  further  research 
illustrating  the  necking  of  thin  sheets  of  aluminum  around  equibiaxially-expanded  holes 
using  a  hydraulic  bulger.  Nazeer  et  al.  (2000  [6])  using  a  conical  tool,  and  Simonsen  et  al. 


(2000  [5])  using  a  spherical  indenter,  analyzed  the  material  mechanics  of  ductile  metal 
sheets. 


Figure  4:  Dishing,  Disking  and  Petalling  of  Plate 
under  (L)  Explosive  Loading  (from  Wierzbicki  [3]); 

(R)  Lateral  Indention  by  a  Sphere  (from  Simonsen  et 
al.  [5]) 

Although  these  studies  were  primarily  concerned  with  relating  ballistic  penetration 
using  quasi-static  methodology,  they  had  a  strong  physical  correlation  with  the  behavior  of 
thin  sheets  subjected  to  dynamic  blast  loading,  see  Figure  4.  Additionally,  after  examining 
the  plate  cutting  behavior  of  vessel  groundings,  Wierzbicki  (et  al.  1993  [7]  and  1999  [3]) 
proposed  that  the  kinematics  of  the  thin  plate  cutting  process,  as  seen  in  Figure  5,  was 
comparable  to  those  of  both  ballistic  penetration  and  explosive  petalling.  To  explore  the 
extent  of  both  of  these  physical  correlations  Woertz  (2002  [4])  studied  the  deformation  of 
clamped  steel  plates  in  two  phases: 

1.  Using  a  spherical  indenter  to  model  early  phase  dishing,  and  subsequendy 
disking. 

2.  Using  an  oblique  conical  punch  to  model  late  phase  radial  crack  propagation 
and  petalling. 


Figure  5:  Similarity  in  the  Kinematics  of  Wedge  Cutting 
(Left)  and  Petalling  (Right)  (from  Wierzbicki  [3]). 

This  approach  was  largely  successful  in  investigating  the  first  phase,  dishing  and 
disking,  but  met  some  difficulty  in  the  second.  In  addition  to  the  physical  limitations  of  the 
equipment  used  to  induce  radial  cracking,  Woertz  also  found  the  frictional  interaction 
between  the  sample  and  the  conical  punch  to  be  problematic  in  analyzing  the  force- 
displacement  history. 


The  effects  of  friction  in  the  wedge  cutting  model  severely  hamper  its  utility  in  the 
friction-free  petalling  phase  of  ballistic  and  blast  failures.  Woertz  assumed  only  two 
components  of  work-energy  dissipation  in  the  petalling  of  thin  metal  sheets,  bending  work 
and  membrane  energy.  Thomas  (1992  [8])  estimated  that  in  addition  to  bending  and 
membrane  work,  friction  accounts  for  as  much  as  40  percent  of  the  work-energy  dissipated 
in  the  mechanics  of  plate  cutting.  Zheng  et  al.  (1996  [11])  characterized  the  frictional  force 
on  a  wedge  in  the  steady-state  cutting  of  a  plate  as  machining  friction,  near  the  tip,  and 
sliding  friction,  along  the  sides  of  the  wedge.  Attempts  by  Lu  et  al.  (1990  [10])  were  made 
to  quantify  this  frictional  component  in  the  cutting  process  by  measuring  the 
disengagement  force  of  the  cutting  wedge.  Yet  no  reliable  method  has  been  developed  to 
accurately  quantify  the  contribution  of  friction  to  the  process  of  wedge  cutting,  and  by 
extension  quasi-statically  model  petalling  and  crack  propagation. 


An  alternative  approach  to  quasi-statically  modeling  crack  propagation  and  petalling 
may  be  to  use  a  variation  of  the  trousers  test  of  tearing  ductile  metal  sheets.  Yu  et  al.  (1988 
[17])  analyzed  the  energy  dissipated  in  bending  and  tearing  thin  aluminum  alloy  sheets  along 


pre-machined  grooves,  using  two  counter-rotating  cylinders,  see  Figure  6.  This  method 
preserved  the  key  elements  of  petalling  kinematics,  including  bending  work  and  membrane 
energy,  but  removed  the  added  effects  of  faction  previously  encountered,  Figure  7. 
However,  by  pre-machining  grooves,  to  guide  the  propagation  of  the  tearing  fracture,  the 
#  material  properties  of  the  sample  were  altered,  affecting  the  results.  Lu  et  al.  (1994  [33]) 

avoided  this  pre-machining  by  fashioning  the  sample  of  thin  metal  plate  into  a  box  column 
and  allowing  the  tearing  fracture  to  propagate  along  the  comers.  This  approach  also 
preserved  the  kinematics  of  petalling,  but  the  geometric  discontinuities  of  the  sample  at  the 
sharp  bends  of  the  comers  may  have  likewise  affected  the  results. 


Figure  6:  Counter-Rotating  Cylinder  Trousers  Test 
(from  Yu  et  al.  [17]) 


Figure  7:  Cylindrical  Roller  Geometry  of  Petalling 
(from  Wierzbicki  [3]) 


A  possible  solution  to  quasi-statically  modeling  the  propagation  of  cracks  and  petalling 
of  thin  plates  builds  upon  the  work  of  Lu  et  al.  Through  the  use  of  a  similarly  configured 
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testing  apparatus  and  test  samples,  with  specific  connection  tab  details,  a  more  accurate 
analysis  of  crack  propagation  and  plate  petalling  may  be  made  that  incorporates  plate 
bending  energy,  and  membrane  energy  but  avoids  the  inclusion  of  frictional,  machining 
and  bending  effects.  This  research  develops  a  detailed  apparatus  design  and  method  to 
conduct  this  analysis  and  compares  testing  results  to  analytically  derived  expected  values. 


TESTING  METHOD  AND  APPARATUS  DESIGN 

DEVELOPMENT 

The  development  of  the  modified  trousers  test  apparatus  and  method,  for  use  as  a 
quasi-static  model  for  crack  propagation  and  petalling  was  conducted  in  three  phases.  The 
first  phase  was  an  analytically  based  investigation  of  crack  propagation  and  petalling  with 
the  purpose  of  defining  gross  load-displacement  requirements  of  a  detailed  testing 
apparatus  design.  The  second  phase  was  a  qualitative  investigation  of  sample  material 
preparation  and  test  apparatus  geometry  in  pursuit  of  an  understanding  of  the 
characteristics,  and  possibly  control  of  fracture  propagation.  The  final  phase  was  to 
develop  a  detailed  testing  apparatus  design  to  be  used  in  future  studies  to  validate  the  quasi¬ 
static  modeling  method  by  testing  samples  of  thin  mild  steel. 

Analytical  Investigation 

Preliminary,  order  of  magnitude,  approximate  analysis  is  included  in  the  Analytical 
Investigation  section  of  this  paper,  below. 

Qualitative  Investigation 

Sample  Preparation 

The  point  of  departure  from  previous  trousers  test  studies  of  this  work  was  the 
specific  geometry  of  the  sample.  Previous  trousers  test  samples  used  flat  metal  plates,  with 
pre-cut,  rectangular  tabs,  tom  in  the  fashion  of  Figure  8.  The  purpose  of  these  tests  was  to 
investigate  the  energy  dissipation  of  tearing  fractures,  not  in  relation  to  cracking  and 
petalling.  Hence,  the  opposite,  “reverse  curvature”  of  every-other  sample  section  was  not 
of  kinematic  concern.  However,  to  relate  this  type  of  tearing  to  crack  propagation  and 
petalling,  including  fracture  and  bending  energy,  it  was  important  to  isolate  the  curvature  to 
a  single  portion  of  the  sample  material. 


Figure  8:  Trousers  Test  Sample  with  Pre-cut 
Rectangular  Tabs  and  Machined  Grooves  (from  Yu 
etal.  [17]). 

To  achieve  this  type  of  tearing  geometry,  the  thin  sample  plates  were  bent  into  box 
columns,  in  the  fashion  of  the  samples  of  Lu  et  al.  However,  while  Lu  attached  the  entire 
box  edge  to  one  of  four  cylindrical  rollers,  pre-cut  tabs,  located  centrally  on  two  opposing 
faces  of  the  column  edge  attached  the  samples  tested  herein.  Figure  9.  This  approach 
isolated  the  bending  and  curvature  induced  by  the  rollers  to  a  flap  of  material  out  of  the 
center  of  two  opposing  faces  of  the  box,  while  maintaining  an  un-curved  geometry  for  the 
remainder  of  the  sample,  better  approximating  the  kinematics  of  petalling. 


Figure  9:  Rough  Sample  Geometry:  (L)  Pre-Cut  Tabs 
Centrally  Located  on  Opposing  Faces,  (R)  Tabs 
Attached  to  Testing  Apparatus. 

In  the  tests  of  Yu  et  al.  (1988  [17])  the  propagation  of  the  tearing  fracture  was 
controlled  with  the  machining  of  grooves.  In  previous  box  column  shaped  samples  the 
bent  comers  of  the  box  column  controlled  the  crack  propagation,  Lu  et  al.  (1994  [33]). 


With  the  samples  of  this  study  there  was  to  be  no  machined  or  geometric  preparations  of 
the  sample  to  govern  the  propagation  of  the  tearing  fracture.  As  a  result,  the  first  phase  of 
this  investigation  was  concerned  with  establishing  the  geometry  of  the  connecting  tabs  on 
the  box  column  sample,  to  best  achieve  data  collection  in  the  third  phase.  The  aim  of  this 
first  phase  of  testing  was  a  qualitative  understanding  of  the  type  of  fracture  propagation  to 
expect  during  further  phases  of  investigation,  to  see  if  the  geometry  of  the  connection  tab 
influences  the  convergence  or  divergence  of  the  fracture  lines. 

Apparatus  Configuration 

In  further  development  of  methods  of  controlling  the  propagation  of  fracture  through 
the  sample  material,  the  effect  of  apparatus  geometry  to  control  the  line  of  fracture 
propagation  through  the  sample  material  was  tested.  Previous  trousers  tests  used 
cylindrical  metal  rollers,  with  smooth  and  parallel  surfaces,  Figure  6.  The  purpose  of  this 
phase  of  testing  was  to  investigate  the  effect  of  altering  the  shape  of  the  surface  of  these 
cylindrical  rollers  and  their  position  relative  to  the  sample  material  to  induce  parallel  lines  of 
fracture  propagation  in  the  sample. 

Apparatus  Design 

The  final  phase  of  this  investigation  consisted  of  utilizing  the  results  of  the  previous 
two  phases  in  the  design  of  a  modified  trousers  testing  apparatus.  This  apparatus  combines 
the  analysis  of  the  sample  pre-cut  and  cylindrical  roller  geometry  to  govern  the  propagation 
of  fractures  through  the  sample  material.  This  apparatus  can  then  be  used  to  conduct  a 
series  of  modified  trousers  tests  on  thin  mild  steel  plate  (up  to  h=lmm)  to  model  the 
petalling  deformation  caused  by  close  proximity  explosions,  and  comparing  the  detailed 
force-displacement  data  collected  and  computed  specific  work  of  fracture  to  analytical 
predictions. 


ANALYTICAL  INVESTIGATION 

The  general  theory  used  in  this  section  was  first  derived  by  Wierzbicki  (1999  [3])  and 
simplified  by  Woertz  (2002  [4]).  They  asserted  that  the  total  work  dissipated  in  cracking 
and  petalling  is  due  to  the  propagation  of  the  radial  cracks,  mechanical  bending  of  the 
petals  and  membrane  deformation.  The  bending  analysis  was  developed  from  mechanical 
relations,  and  the  membrane  deformation  derivation  is  an  extension  of  the  derivations  of 
Wierzbicki  et  al.  (1993  [7]).  A  full  derivation  of  force-displacement  relations  is  included  in 
Appendix  A. 

General  Petalling 

Begin  from  very  general  petalling  geometry  of  n  cracks  propagating  from  a  single 
point,  dividing  a  thin  plate  into  n  symmetric  petals,  Figure  10. 


Figure  10:  Theoretical  Petalling  Geometry 


The  central  angle  of  each  petal  is  defined  as  20,  such  that 
0=* 

n  (1) 
And  each  petal  can  be  described  as  a  triangle  OAB.  The  instantaneous  length  of  the 

crack,  X,  is  related  to  the  total  petal  length,  A. 


Bending  Energy 

As  the  petal  grows  in  size,  and  the  radial  cracks  propagate  through  the  material,  the 
hinge  line  AB  moves  through  the  material,  leaving  the  curled  petal  behind.  This  kinematic 
boundary  condition  imposes  a  relation  between  the  propagation  speed  of  this  hinge  line, 
dA/dt,  the  instantaneous  rate  of  rotation  of  the  petal  at  the  hinge  line,  d(|>/dt,  and  the 
instantaneous  petal  radius  of  curvature,  p. 


Wierzbicki  ultimately  derived  an  expression  for  p.  In  this  study,  the  characteristics  of 
the  fixed  cylinders  of  the  testing  apparatus  dictate  that  the  instantaneous  radius  of  curvature 
of  the  petals  is  known  and  constant. 


Continuing  the  assumption  of  a  rigid,  perfectly  plastic  material,  with  an  average  flow 
stress  of  O,  the  fully  plastic  bending  moment  per  unit  length  of  the  flat  metal  sheet,  using 


the  Tresca  yield  criteria,  is: 


(4) 


where  h  is  the  plate  thickness.  Although  Wierzbicki  and  Woertz  continued  to  state  that  the 
curved,  dished  surface  of  the  thin  plate  would  stiffen,  and  amplify  the  plastic  bending 
moment  by  the  amplification  factor  r|,  the  thin  plate  of  this  study  remained  flat  and  un¬ 
dished.  Hence,  T| = 1 . 


The  rate  of  bending  work  of  one  petal  is  expressed  as: 

iwb  =  2.m0.lab.1* 


where  LAB— 2Atan9.  Substituting  Equation  (3)  into  Equation  (5)  yields: 


To  apply  this  to  the  quasi-static  model  of  petalling  used  in  this  study,  the  changes  in 
work  dissipated  over  short  increments  of  time  and  small  increments  of  displacement  can  be 
obtained  by  integrating  in  time: 


and  Equation  (6)  becomes: 


4-M  •  A 

Wb  = - tan(e) 

P 


Tearing  Energy 

For  perfectly  britde  materials,  the  crack  width  between  adjacent  petals  can  be 
expressed  as  a  function  of  the  distance  from  the  point  of  intersection  of  two  adjacent  hinge 
lines. 

,3 

8(x)  =  — — -sin(0)-cos(0) 

3  2 

P  < 
where  8  is  the  local  crack  width,  x  is  very  near  the  crack  tip,  and  p  is  constant  In  real, 
ductile  material,  the  crack  tip  does  not  coincide  with  the  intersection  of  the  hinge  lines,  but 
where  local  strain  reaches  the  crack  tip  opening  displacement  parameter  (CTOD) 
(Wierzbicki  et  al.  1993  [7]),  5t.  The  length  of  the  plastic  zone  near  the  crack  tip  can  be 
found  using  CTOD  and  Equation  (9): 


a  o  Q  _  1 

Xp  =  1.44p  -8t  -sin(0)  •cos  (e)~  ‘  (l 

Leading  to  the  calculation  of  the  rate  of  membrane  energy  dissipation  in  the  plastic 
zone,  near  the  crack  tip: 


2  U  d  A 

7CTohy-A 
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Using  Equations  (4)  and  (10)  in  Equation  (11): 


-W  = 
.  m 


3.84M  *6t  -p  -sin(0)  —  A 

1 _ dt 

h-cos(0) 


To  apply  this  to  the  same  quasi-static  model  of  petalling  used  in  this  study,  the  changes 
in  work  dissipated  over  short  increments  of  time  and  small  increments  of  displacement 
were  again  estimated,  from  Equation  (7),  and  Equation  (12)  becomes: 


W  = 
m 


3.84M0-8t  p  -A-sin(e) 
h'cos(0) 


Total  Energy 

Adding  Equations  (8)  and  (13)  to  get  the  total  energy: 


wt=  wb  + wm 


4-Mq-  A2-tan(0)  3.84MQ- A-5t  3 -p  3  *sin(0)  3 

Wt=  -  + - - 

p  h-cos(0j  q 

To  apply  this  to  the  quasi-static  model  of  petalling  used  in  this  study,  the  changes  in 
work  dissipated  over  short  increments  of  time  and  small  increments  of  displacement  were 
again  estimated  and  Equation  (15)  becomes: 

F  =  —  W4 


A  force-displacement  trace  for  this  expression  was  generated  for  comparison  between  this 
analysis,  and  wedge  cutting  analysis,  and  is  included  in  Appendix  B.  A  general  example  of 
the  generated  force-displacement  curves  is  computed  in  Appendix  A,  and  included  in 
Figure  11. 
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Force  v.  Petal  Length 


X 

Crack  Length  (m) 

Figure  11:  Approximate  Theoretical  Load- 
Displacement  Curve  for  Petalled  Plate 


Expected  Sample  Energies 

For  the  samples  tested  in  this  study,  the  geometry  does  not  follow  the  general  petalling 
geometry.  Without  the  geometric  or  machined  details  of  previous  trousers  tests  the 
propagation  of  the  fractures  follows  a  path  similar  to  those  described  in  Simonsen  et  al. 
(1997  [12])  for  the  concertina  tearing  mode  of  plate  failure.  That  is,  the  fracture 
propagation  lines  will  not  follow  the  angular  petal  lines,  but  will  either  become  convergent 
or  divergent,  Figure  12.  -■  ' 
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Figure  12:  (L)  Converging  Fracture  and  (R) 
Diverging  Fracture  Geometries. 


To  that  end,  the  sample  and  testing  apparatus  geometry  were  set  to  induce  nearly 
parallel  fractures.  Hence,  the  general  fracture  geometry  is  no  longer  the  triangular  petals 
previously  discussed,  but  becomes  that  of  Figure  13. 


Figure  13:  Sample  Petalling  Geometry 


The  lines  of  propagation  of  the  non-ideal  petal,  or  tab,  are  idealized  as  parallel  and 
each  petal  can  be  described  as  a  rectangular  tab.  The  instantaneous  length  of  the  crack,  X, 
can  be  defined  as  a  function  of  the  total  petal/tab  length,  A. 


X-  A  -  C 


(17) 


Where  C  is  the  pre-cut  length. 


Bending  Energy 

Using  the  same  assumption  of  a  rigid,  perfecdy  plastic  material,  with  an  average  flow 

stress  of  Cc,  then  the  fully  plastic  bending  moment  per  unit  length  of  the  flat  metal  sheet, 

27 


m 


using  the  Tresca  yield  criteria,  was  previously  derived,  as  Equation  (4).  The  rate  of  bending 
work  of  one  petal  is  expressed  in  Equation  (5)  with  L^^b,  or: 


— Wh  =  2M_b  —  <(> 

1  dt  D  0  dt 

Substituting  Equation  (3)  into  Equation  (18)  yields: 


(18) 


=  2-M0.b- 


To  apply  this  to  the  quasi-static  model  of  petalling  used  in  this  study,  the  changes  in 
work  dissipated  over  short  increments  of  time  and  small  increments  of  displacement  can  be 
obtained  by  integrating  Equation  (19)  in  time,  to  become: 


2-M  -b-A 

wb= - 


(20) 


♦ 


# 


Tearing  Energy 

Continuing  the  assumptions  of  the  previous  analysis,  for  perfectly  brittle  materials,  the 
membrane  energy  rate  of  dissipation  of  on  petal  remains  unchanged  from  Equation  (12). 
The  total  membrane  energy  remains  unchanged  from  Equation  (13),  with  0=60  degrees. 


Total  Energy 

Adding  Equations  (20)  and  (13)  to  get  the  total  energy: 

J_  _2 

lQ-b A  3.84M0  A-5t3  p3  sin(e)  3 
-  +  - 


Wt  = 


h-cos  (©) 


(21) 


To  apply  this  to  the  quasi-static  model  of  petalling  used  in  this  study,  the  changes  in 
work  dissipated  over  short  increments  of  time  and  small  increments  of  displacement  were 
again  calculated  and  Equation  (21)  becomes: 


(22) 
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Force  (N) 


# 


* 


A  general  example  of  the  generated  force-displacement  curves  is  computed  in 
Appendix  C,  and  included  in  Figure  14.  Force-displacement  curves  corresponding  to  each 
sample  tested  are  computed  and  included  in  Appendix  H. 


Figure  14:  Approximate  Theoretical  Load- 
Displacement  Curve  for  Tabbed/Petalled  Sample 
Plate 


# 


# 
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QUALITATIVE  INVESTIGATION 


Sample  Preparation 

Method 

r 

To  accomplish  the  second  phase  of  this  study,  qualitative  investigations  of  various  tab 
geometries  were  carried  out  on  thin  gauge  aluminum  sheet,  h=0.1 117mm.  Five  tab 
geometries  were  fabricated  onto  the  edges  of  flat  samples,  Figure  15,  clamped  on  all  four 
sides.  The  samples  were  subjected  to  tearing  fractures  using  a  rolling  cylinder  of  radius 
p=  1.5cm,  and  the  behavior  of  the  fracture  propagation  was  noted  and  photographed, 
Appendix  D. 

(a)  (b)  (c) 


60* 

^ —  "~''y 

V 

so* 

EJ3 

^ . . . '"'V 

\  / 

<50’ 

;  \  / 

(d)  (e) 


Figure  15:  Qualitative  Tab  Sample  Geometry  (a)  Six 
Petal  Configuration,  (b)  Four  Petal,  (c)  Trousers 
Configuration,  (d)  Six  Petal  Wide  Tab  Configuration, 

(e)  Four  Petal  Wide  Tab  Configuration. 

The  first  sample  preparation  was  fabricated  with  pre-cut  notches  inclined  at  60  degrees 
from  the  free  edge,  forming  a  60-degree,  triangular  tab.  This  configuration  was  included  to 
reproduce  the  geometry  of  a  six-petal  blast  hole.  The  second  sample  was  fabricated  with 
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pre-cut  notches  inclined  at  45  degrees  from  the  free  edge,  forming  a  90-degree,  triangular 
tab.  This  configuration  was  to  reproduce  the  geometry  of  a  four-petal  blast  hole.  The  third 
sample  configuration  was  fabricated  with  two  parallel,  pre-cut  notches,  forming  a 
rectangular  tab.  This  configuration  was  used  for  comparison  to  standard  trousers  test 
geometries.  The  fourth  and  fifth  configurations  were  to  reproduce  the  six  and  four  petal 
geometry,  respectively,  with  wider  tabs  to  possibly  accommodate  fracture  propagation. 

Results 

Figure  16  illustrates  three  cases  of  the  fracture  geometry  encountered  in  the  first  phase 
of  testing.  The  complete  results  of  this  testing  phase  are  located  in  Appendix  D,  The 
sample  geometries  shown  in  the  figure  are  the  six  petal  pre-cut,  the  six  petal  wide  tab,  and 
the  parallel  pre-cut  tab  arrangement.  From  these  representative  cases  it  is  seen  that  the 
cracks  followed  neither  the  line  of  the  angled  pre-cuts  nor  ran  parallel  through  the 
aluminum  sheet. 


Figure  16:  Phase  One  Results  (L  to  R)  n=6  Pre-cut 
Tab,  n=6  Wide  Pre-cut  Tab  and  Parallel  Pre-cut 
Tab. 


As  expected,  all  samples  tested  exhibited  converging  fracture  lines,  independent  of  pre¬ 
cut  tab  geometry.  Further,  the  cracks  of  all  samples  converged  at  a  relatively  shallow  angle, 
ranging  from  7  to  10  degrees,  and  remained  fairly  straight.  Complete  results  are  found  in 
Figure  17. 


Figure  17:  Comprehensive  P' 

aase  One  Results 

Sample 

Number 

Pre-Cut 

Geometry 

Maximum 
Tab  Width 
(mm) 

Average  Fracture 
Convergence 
Angle  (deg) 

Total  Effective 
Fracture 
Length  (mm) 

i 

HUSH 

30 

6.8 

80 

2 

50 

9.25 

80 

3 

HESS 

11.7 

10.15 

41 

4 

50 

7.8 

80 

5 

Parallel 

50 

8.25 

80 

# 


m 
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Discussion 

From  the  results  of  this  first  phase  of  investigation  it  was  seen  that  the  line  of 
propagation  of  the  fractures  induced  by  a  rolling  cylinder  were  independent  of  the  pre-cut 
tab  geometry.  The  angled  pre-cuts  experienced  converging  fracture  lines  of  very  similar 
convergence  angles  as  the  parallel  pre-cuts.  Further,  all  of  the  samples  exhibited  fairly 
constant  convergence  angles,  resulting  in  straight  fracture  lines. 

It  was  also  observed  in  this  phase  of  investigation  that  although  the  thin  aluminum 
sheet  was  tightly  clamped  as  the  rolling  cylinder  progressed,  the  sample  was  stretched  and 
became  raised,  or  bowed,  in  the  region  of  the  rolled  tab.  It  may  be  this  bowing  curvature 
and  stretching  of  the  material  that  induced  the  converging  fracture  geometry.  It  was  this 
observation  that  provided  motivation  to  conduct  the  second  phase  of  investigation. 

As  a  result  of  this  first  phase,  it  is  asserted  that  the  pre-cuts  in  the  boxed  material 
samples  should  be  fabricated  to  ease  connection  of  the  sample  to  the  testing  apparatus,  and 
maximize  the  overall  fracture  length.  Both  objectives  may  be  achieved  by  widely  spacing 
the  pre-cuts  on  the  face  of  the  sample.  The  wider  tab  allows  for  a  more  secure  connection 
between  the  sample  material  and  the  surface  of  the  apparatus.  The  wider  tab  also  allows 
for  a  longer  fracture  length  before  the  fracture  lines  converge  upon  each  other. 

Additionally,  the  fairly  constant  angle  of  convergence  encountered  in  this  phase  of 
testing  suggests  that  analytical  approximations  of  the  force-displacement  relations  for  each 
sample  may  be  improved  to  account  for  this  non-ideal  fracture  line  geometry.  The  exact 
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convergence  angle  of  each  sample  may  be  measured  to  impose  this  correction,  or  for  very 
Shallowly  converging  cracks  the  fracture  line  may  be  approximated  as  parallel  for  analysis. 

Apparatus  Configuration 

Method 

To  accomplish  the  second  phase  of  testing,  qualitative  investigations  of  various  rolling 
cylinder  geometries  were  carried  out  on  thin  gauge  aluminum  sheet,  h=0.1117mm,  with 
wide  tab,  60  degree  pre-cuts  (Figure  15d).  The  samples  were  subjected  to  tearing  fractures 
using  three  rolling  cylinder  face  geometries,  Figure  18,  as  in  phase  one  of  testing. 
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Figure  1 8:  Qualitative  Rolling  Cylinder  Geometry;  (a) 

Parallel  Cylinder,  (b)  Conically  Tapered  Cylinder,  (c) 

Spherically  Tapered  Cylinder. 

The  first  cylinder  tested  was  a  simple,  parallel  roller  of  p=15mm.  This  configuration 
was  included  to  reproduce  and  compare  the  results  encountered  in  the  first  phase  of  testing. 
The  second  cylinder  was  fabricated  with  a  conically  tapering  radius,  pmax=20mm.  The  third 
cylinder  configuration  was  fabricated  with  a  spherically  tapering  radius,  pmax=20mm, 
PsPh«e=65mm. 

The  two  tapered  cylinders  were  connected  to  the  tabs  of  the  thin  aluminum  in  two 
configurations,  Figure  19,  flush  to  the  point  of  maximum  radius  and  recessed  to  the  15mm 
uniform  radius. 


# 


33 


# 


Figure  19:  Phase  Two  Connection  Geometry  (a) 

Flush,  (b)  Recessed 

The  samples  were  subjected  to  tearing  fractures  using  the  three  rolling  cylinders,  and 
the  two  connection  geometries  and  the  behavior  of  the  fracture  propagation  was  noted  and 
photographed,  Appendix  F 

Results 

Figure  20  illustrates  three  cases  of  the  fracture  geometry  encountered  in  the  second 
phase  of  testing.  The  complete  results  of  this  testing  phase  are  located  in  Appendix  F.  The 
sample  geometry  used  in  this  phase,  and  shown  in  the  figure,  was  the  six  petal  pre-cut  wide 
tab  arrangement.  The  apparatus  geometries  illustrated  in  the  figure  are  of  the  parallel  and 
conically  tapered  cylinders,  in  the  flush  and  recess  mounted  configurations.  From  these 
representative  cases  it  is  seen  that  the  path  of  crack  propagation  was  influenced  by  the 
geometry  of  the  sample  rolling  apparatus.  a 


Figure  20:  Phase  Two  Results  (L  to  R)  Parallel  Face 
Cylinder,  Flush  Mounted  Conically  Tapered  Cylinder 
and  Recess  Mounted  Conically  Tapered  Cylinder. 
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As  seen  in  the  complete  Phase  Two  data,  the  geometry  of  the  cylinder  face  alone  did 
not  have  a  profound  effect  upon  the  convergence  of  the  fracture  lines  through  the  thin 
aluminum  sample.  That  is,  the  flush  mounted  conically  tapered  cylinder  had  fracture  lines 
converging  at  a  rate  not  dissimilar  to  those  of  spherically  tapered  cylinder  and  the  parallel, 
simple  cylinder.  Of  greater  effect  upon  the  fracture  propagation  was  the  detail  of 
connection  between  the  sample  and  the  roller.  Specifically,  recessing  the  conically  or 
spherically  tapering  segment  of  the  cylinder  significantly  reduced  the  angle  of  convergence 
of  the  ffactute  lines  in  the  samples.  Complete  results  are  found  in  Figure  21 . 


Figure  21:  Comprehensive  Phase  Two  Results 


Sample 

Number 

Cylinder 

Geometry 

Connection 

Geometry 

Average  Fracture 
Convergence 
Angle  (deg) 

1 

Parallel 

Flush 

Mounted 

7.8 

2 

Conically 

Tapering 

Flush 

Mounted 

5.65 

3 

Spherically 

Tapering 

Flush 

Mounted 

4.55 

4 

Conically 

Tapering 

Recessed 

2.1 

5 

Spherically 

Tapering 

Recessed 

2.5 

Discussion 

From  the  results  of  this  second  phase  of  investigation  it  was  seen  that  the  line  of 
propagation  of  the  fractures  induced  by  a  rolling  cylinder  were  influenced  by  the  geometry 
of  the  face  of  the  cylinder.  Cylinders  with  regions  of  convex  tapered  radii  induced 
shallower  angles  of  fracture  convergence  than  simple,  parallel-faced  cylinders.  Further,  it 
was  found  that  the  method  of  connection  between  the  thin  aluminum  sample  and  the 
convex  tapered  cylinder  also  influenced  the  angle  of  convergence.  By  recessing  the  region 
of  convexity  into  the  sample,  the  result  of  connecting  the  sample  material  to  the  non- 
convex  parallel  region  of  the  cylindrical  roller,  the  fracture  lines  could  be  made  nearly 
parallel. 


As  a  result,  it  is  asserted  that  the  propagation  of  the  fracture  lines  through  the  sample 
material  in  the  third  phase  of  this  investigation  may  be  controlled  using  methods  other  than 
physically  altering  the  sample.  Through  the  addition  of  a  convex  region  to  the  face  of  the 
rolling  cylinders,  and  the  recessed  attachment  of  the  sample  material  to  the  cylinders,  nearly 
parallel  fracture  lines  can  be  induced. 


APPARATUS  DESIGN 


Using  the  results  of  the  second  phase  of  testing,  box  column  samples  as  seen  in  Figure 
22  are  to  be  constructed  for  use  on  the  test  apparatus.  Complete  test  sample  specifications 
are  included  in  Appendix  E.  The  samples  are  to  be  inserted  into  the  testing  apparatus  as 
illustrated  in  Figure  23.  Complete  testing  apparatus  specifications  are  included  in  Appendix 
G.  The  two  rollers  are  then  to  be  driven  simultaneously  by  pulling  up  the  four  attached 
wire  ropes.  This  motion  causes  bending  of  the  two  pre-cut  tabs  onto  the  rollers,  and  at  the 
same  time  propagates  tearing  along  the  two  opposite  sample  faces. 


Figure  22:  Box  Column  Sample  Geometry 


Figure  23:  Schematic  of  Experimental  Setup 


CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions: 

Through  the  preliminary  theoretical  analysis  it  was  found  that  comparison  between 
wedge  cutting  and  petalling  kinematics,  based  on  similar  geometry,  is  well  grounded.  As 
seen  in  Appendix  B,  force  displacement  relationships  in  each  case  are  nearly  identical  when 
the  effect  of  wedge  cutting  friction  is  neglected.  However,  in  actual  wedge  cutting 
processes,  the  effect  of  friction  is  great,  and  negates  the  utility  of  comparison  between  the 
two  phenomena. 

As  a  superior  method  of  recreating  the  same  petalling  process  using  a  quasi-static 
approach,  the  modified  trousers  test  of  this  study  proved  useful.  It  was  initially  found  that 
petalling-like  fractures  in  thin  metal  samples  could  be  reliably  produced  using  a  rolling 
cylinder,  which  generated  converging  lines  of  fracture.  This  convergence  was  postulated  to 
be  the  result  of  the  ductile  characteristics  of  the  thin  samples  giving  rise  to  mixed  mode,  in¬ 
plane  and  out-of-plane,  tearing.  Hence,  analyzing  force-displacement  data  using  an 
idealized,  triangular  petal  was  revealed  to  be  inadequate.  The  plastic  hinge  line,  propagating 
away  from  the  tip  of  the  petal,  decreases  in  length  in  a  converging  geometry,  as  opposed  to 
expanding  in  length  in  the  idealized  model. 

It  proved  more  accurate  to  model  the  petalling  fracture  propagation  as  a  rectangular 
tab,  with  parallel  lines  of  fracture.  This  compelled  the  development  of  a  method  to  offset 
the  convergence  of  fracture  lines,  and  produce  parallel  fractures.  This  was  to  be  achieved 
without  altering  the  sample  geometry  or  material  properties,  as  inadvertently  done  in  past 
trousers  tests  (Lu  et  al.  1994  [33]  and  Yu  et  al.  1988  [17]).  It  was  shown  that  altering  the 
geometry  of  the  sample  pre-cuts  had  little  effect  on  the  lines  of  fracture  through  the 
sample.  Better  results  were  achieved  in  controlling  the  fracture  convergence  by  changing 
the  shape  of  the  cylindrical  roller.  Through  the  addition  of  a  raised  portion  to  the  cylinder 
face,  the  angles  of  convergence  in  the  sample  were  modestly  reduced.  Altering  the 
connection  geometry  of  the  raised  portion  of  the  cylinder  to  the  sample  material  proved 
most  successful  in  controlling  the  angle  of  fracture  convergence.  Through  the  combination 
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of  the  raised  cylinder  face  and  modified  connection  detail,  the  lines  of  fracture  were  made 
nearly  parallel. 

With  this  knowledge,  the  modified  trousers  testing  apparatus  was  constructed,  as 
detailed  in  Appendix  G,  and  sample  specimens  were  fabricated,  as  detailed  in  Appendix  E. 

Recommendations 

Most  necessary  in  future  work  is  an  extensive  study  using  the  newly  constructed  testing 
apparatus  and  modified  trousers  test  method  to  thoroughly  validate  this  quasi-static 
approximation  of  the  fracture  and  petalling  phenomena.  Further  investigations  should 
include  a  range  of  sample  thickness,  up  to  a  maximum  thickness  of  1mm  (for  mild  steel). 
Additionally,  through  testing  multiple  materials  a  greater  understanding  of  the  petalling 
process  could  be  achieved. 

In  future  studies  using  the  modified  trousers  test  apparatus  an  analysis  of  the  strain 
field,  near  the  crack  tip  should  be  conducted.  In  a  method  similar  to  that  used  by  Woertz 
(2002  [4]),  fine  grid  markings  could  be  made  on  the  sample  faces  to  compute  instantaneous 
local  strains.  These  strains  could  be  used  as  a  measure  of  material  stretching,  three 
dimensional  petal  displacements,  and  bending  work  dissipated  in  achieving  the  final 
deformed  geometry.  From  such  an  analysis,  a  more  defined  understanding  of  total  work 
dissipated  could  be  achieved,  aiding  in  an  understanding  of  the  methods  of  energy 
dissipation  in  the  initial  phase  of  dynamic,  explosive  events. 

\ 

In  conjunction  with  an  investigation  of  the  strain  field,  perhaps  a  superior  method  of 
computing  the  displacement  of  the  sample  material  could  be  developed.  In  stead  of  wholly 
relying  on  the  position  of  the  testing  machine  cross-head,  and  from  there  calculating  the 
relative  motion  of  the  sample,  to  determine  strain,  measuring  displacements  directly  from 
the  sample  may  prove  more  useful  for  the  analysis  of  the  crack  tip  strain  field.  Measuring 
sample  displacements  directly  could  be  effected  indirectly,  by  measuring  the  angle  of 
rotation  of  each  roller,  or  direcdy,  by  measuring  a  pre-determined  point  on  the  sample  to  a 
fixed  point  in  space. 


39 


Perhaps  most  significantly,  results  from  future  tests  using  the  modified  trousers  test 
apparatus  should  be  compared  to  results  obtained  from  existing  numerical  models.  Such 
models  could  be  constructed  in  ABAQUS  or  LS-DYNA  to  investigate  the  mode  and 
location  of  fracture,  and  approximate  the  deformed  sample  shape  as  the  result  of  fracture. 
Further,  such  numerical  results  would  serve  to  improve  upon  the  design  of  the  apparatus. 
Specifically,  through  an  understanding  of  the  details  of  fracture  the  shape  of  the  cylindrical 
rollers  could  be  improved  to  further  control  fracture  convergence. 
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APPENDIX  A:  PETALLING  FORCE-DISPLACEMENT 

APPROXIMATION 

For  a  sample  plate  of  thin,  ductile  metal  with  the  following  characteristics: 

Plate  thickness 


h  :=  .5 -mm 


G0  :=  27210  -Pa 


CTOA  :=  10- deg 


0  :=  30- deg 


AQ  :=  L5-cm 


Average  Flow  Stress 

Crack  tip  opening  angle  (CTOA) 


Corresponding  to  petal  semi-angle 
where  n=6 


Pre-cut  petal  length 
On  the  testing  apparatus  with  the  following  characteristics: 


P0  :=  1.5-cm 


Rolling  cylinder  radius 


p  j  :=  1.25 -cm 


Adot  :=  10 


mm 

min 


a(a.)  :=  ^-cos(0)  +  A 


A.(a)  :=  A — “‘Cos (0)  1 
Pwr 


Resulting  in: 


A  A  P° 

Adot  :==  Adot - 

Pi 

5ct0d  M  :=  2-X-sin(CTOA) 


Po 


5ctod(A)  :=  2’ A — —  *cos(0)  1  sin(CTOA) 
Pwr 


Wire  rope  cylinder  radius 
Cross-Head  vertical  speed 

Total  petal  length  as  a  function  of 
fracture  length 

Fracture  length  as  a  function  of 
cross-head  vertical  displacement 


Petal  length  rate  of  change 

Crack  tip  opening  distance  as  a 
function  of  fracture  length 

CTOD  as  a  function  of  cross-head 


Total  bending  moment  per  petal  per  unit  length 
2 


M0  := 


Vh 


wb(x)  := 


4-M0-(a(a)  -  A0)2  tan(e) 
Po 


displacement 


Mq  =  17- 


N-m 


m 


Total  bending  work  per  petal  as  a 
function  of  fracture  length 


Total  bending  work  per  petal  as  a 


4-M0- 


Wb(A)  := 


V  Pwr  J 


■tan(0) 


Po 

function  of  cross-head  displacement 
And  the  contribution  of  membrane  work  was  expressed  as: 


WmM  :=  M0-(aU)  -  A0)-3.84h  ti8ctodM)  J'(Po)  J'sin(9)  3  -cos(e) 


1  2  -4 

3  /  \  3  .  /„\  3  1 


or: 


r  ~  \ 


I  1  -4 

3  /  \3  .  3  /^y-  1 


Wm(A):^Mo-  A-—  *3.84h  l{5ctod(A)j' 3 .(p0) 3 -sin(e)  3  -cos(e) 


y  Pwry 


Making  the  total  work: 

Wt(A)  :=  Wm(A)  +  Wb(A) 
And  the  total  force: 


f(a)  :=  —  Wt(A) 
dA 


\ 
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APPENDIX  B:  PF.TAT  I  .TNr,  AND  WEDGE  CUTTING 


For  a  sample  plate  of  thin,  ductile  metal  with  the  following  characteristics: 


h  :=  .5’ mu 


Plate  thickness 


or,  :=  272-1 06-Pa  ■  _ 

Average  Flow  Stress 

CTOA  :=  10- deg  Crack  tip  opening  angle  (CTOA) 

And  petalling  geometry: 


0  :=  30- deg 


A(A,)  :=  A.-cos(g) 


Corresponding  to  petal  semi-angle 
where  n=6 

Petal  length  as  a  function  of  crack 
length 


Woertz  (2002  [4])  built  upon  the  derivations  of  Wierzbicki  (1999  [3])  to  derive 
simplified  expressions  for  the  total  work  dissipated  in  the  formation  of  radial  cracks  and 
petals.  For  application  to  the  testing  of  this  work,  the  instantaneous  radius  of  curvature  of 


the  petals  was  made  constant,  to  reflect  the  curvature  induced  by  the  cylinders  of  the  testing 
apparatus: 


P0  :=  ’1.5  cm 

Other  characteristics  of  the  testing  apparatus: 


Adot  :=  10 


mm 

min 


Rolling  cylinder  radius 


Cross  head  speed 


a{x)  :=  a(a) 


Cross  head  vertical  displacement 

Woertz  also  decomposed  the  total  work  into  contributions  of  bending  work  and 
membrane  work.  As  that  petalling  is  a  frictionless  process,  he  included  no  contribution  of 
friction  in  his  simplified  expressions.  The  bending  work  was  expressed  as: 


:= 


<Vh 


M0  =  17 


Nm 


m 


Total  bending  moment  per  petal  per  unit  length 

4M  -(A(x.))2-tan(e) 

■  Wb(x)  := - 

^  0  Total  bending  work  per  petal 

And  the  contribution  of  membrane  work  was  expressed  as: 


5ctod  (A.)  :=  2- A.- sin  (CTO A) 


Crack  tip  opening  distance 


Wm(x)  :=  M0  A(A.)-3.84h  1-(5ctodW)'-(po)3  sin(0)  .  *  -cos^) 


1  2  —  4 

3  /  \3  ,  /A\  3  1 


Making  the  total  work: 

WtW  :=  Wm(A.)  +  Wb(x) 
And  the  total  force: 

fM  :=— Wt(x) 

dA 
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For  a  petal  length  of  70mm  the  force  and  work  are: 

F(70mni)  =  8.004x  103  N  Wt(70-mm)  =  415.384J 

And  traces  of  force  and  work  as  a  function  of  crack  length: 
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X 

Crack  Length  (mm) 

For  analytic  comparison,  Wierzbicki  &  Thomas  (1993  [7])  derive  expressions  to 
produce  the  minimum  cutting  force  of  a  wedge  through  a  thin  plate  with  the  following 
characteristics: 


Non-dimensional  CTOD  parameter 


^  wedge  *  ^eg  Wedge  semi-angle  equal  to  the 

petalling  angle,  corresponding  to 
n=6 

As  the  sum  of  three  components: 

F  =  Fu  +  F  +  F.f 
w  1  b  ^  m  ^  f 

Where: 

Fw  =  Minimum  Cutting  Force  for  One  Fracture 
Fb  =  Flap  Bending  Force  for  One  Fracture 


Fm  =  Membrane  Force  for  One  Fracture 

Ff  =  Friction  Force  for  One  Fracture 

The  underlying  assertion  of  Wierzbicki  &  Thomas  is  that,  with  the  elimination  of 
wedge  friction,  accounted  for  in  their  derivation  but  difficult  to  experimentally  measure,  the 
wedge  cutting  model  can  be  successfully  applied  to  the  petalling  and  cracking  model. 
Hence,  for  purposes  of  comparison  to  Woertz,  the  frictional  component  of  the  crack 
propagation  is  ignored,  and  the  expression  derived  is: 

FwW  •=  1.67-cr0*8mt(A,)  -h  -X  •  sin wedge  )  ’cos  wedge  ) 


Leading  to  the  derivation  of  the  work  dissipated  in  one  fracture  as  a  function  of 
fracture  length: 


:=  Fw((|>)  d<t> 


To  apply  these  expressions  for  use  in  crack  propagation  and  petalling,  it  is  most 
important  to  notice  that  each  petal  consists  of  two  of  these  wedge-like  fractures.  Hence: 

FWtW  :=  2*Fw(a)  Total  Petalling  Force  (Wierzbicki  & 

Thomas)  as  a  function  of  theoretical 
petal  length 

WWT W  :=  2-WtwM  Total  Petalling  Work  (Wierzbicki  & 

Thomas)  as  a  function  of  theoretical 
petal  length 

For  a  theoretical  petal  length  of  70mm  the  force  and  work  are: 


FWT(70-mm)  =  7.733  x  103N 


WWT(70mm)  =  338.297T 


And  traces  of  force  and  work  as  a  function  of  theoretical  petal  length: 


These  results  resemble  data  from  many  wedge  cutting  studies.  The  Force  plot 
indicates  that  with  the  application  of  a  small  force  on  the  wedge  a  small  cut  is  initiated,  and 
that  there  is  no  threshold  force  required  for  the  onset  of  fracture. 

For  material  with  the  same  properties,  and  with  equivalent  fracture  lengths,  it  is  seen 
that  the  wedge  derivations  of  minimum  cutting  force  of  Wierzbicki  et  al.  do  not  compare 
well  with  the  simplified  petalling  derivations  of  Woertz.  With  respect  to  force,  Woertz’s 
petalling  derivations  seem  much  more  likely,  indicating  a  minimum  force  before  petalling  is 
initiated  and  cracks  begin  to  form.  The  wedge  derivations  indicate  that  fracture  occurs 
almost  immediately  upon  application  of  force. 


Wedge  Minimum  Cutting  Force 
•  '  *  •  Radial  Cracking  and  Petalling  Force 


A  comparison  of  the  derivations  of  work  dissipated  in  each  case  proves  to  be  similar. 
It  is  seen  that  the  derivation  of  work  for  the  wedge  cutting  and  petalling  processes  compare 
well,  for  samples  of  similar  characteristics. 


APPENDIX  C:  TABBING/PETALLING  FORCE- 
DISPLACEMENT  APPROXIMATION 


a  sample  plate  of  thin,  ductile  metal  with  the 
h  :=  .5  mm 

ct0  :=  272-106-Pa 


following  characteristics: 
Plate  thickness 

Average  Flow  Stress 


0  Pre-cut  tab / petal  length 

CTOA  :=  10- deg  Crack  tip  opening  angle  (CTOA) 

And  tab/petalling  geometry: 


9  :=  30- deg 


Corresponding  to  petal  semi-angle 


where  n— 6 


b  :=  3  •  cm 


Approximately  constant  tab/petal 


width 


On  the  testing  apparatus  with  the  following  characteristics: 


p  0  :=  1.5  •  cm 


Rolling  cylinder  radius 


p  :=  3 -cm 


Wire  rope  cylinder  radius 


Adot  :=  !<►- 


Cross  head  vertical  speed 


a(a,)  i—  A,  +  A, 


Total  petal  length  as  a  function  of 


fracture  length 


^(a)  :=  A — — 
P  wr 


Fracture  length  as  a  function  of 


Resulting  in: 


cross-head  displacement 


Adot :=  Adof 


Petal  length  rate  of  change 


5ctodU)  :=  2-A.-sin(CTOA) 


Crack  tip  opening  distance  as  a 


function  of  fracture  length 


6ctod(A)  ;=  2‘A — — sin(CTOA) 
Pwr 


CTOD  as  a  function  of  cross-head 


Total  beading  moment  per  petal  per  unit  length 


displacement 


M0  := 


M0  =  17- 
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WbM  := 


2-M0-(a(\)  -  A0)  b 


Total  bending  work  pet  petal  as  a 
function  of  fracture  length 


2-M. 


Wjil  = 


V  pwr  J 


•b 


Total  bending  work  as  a  function  of 


cross-head  displacement 

And  the  contribution  of  membrane  work  was  expressed  as: 

1  2  -4, 

Wm(x)  :=  Mo  (a00  -  A0).3.84h~  3 -(p0) 3 -sin(e)  3  cos(e)" 


\-l 


or: 


<  Po^ 


1  2  -  4 

3  /  \  3  .  /  \  3  1 


Wm(A):=M0-A - 3.84h  -^(a))  *(p0)  -sin(e)  -cos(q) 


\  Pwry 
Making  the  total  work: 

Wt(A)  :=  Wm(A)  +  Wb(A) 
And  the  total  force: 


f(a)  :=  —  Wj(a) 

dA 
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Work  (J)  Force  (N) 


Traces  of  force  and  work  as  a  function  of  cross-head  displacement: 


A 

Cross-head  Displacement  (mm) 


APPENDIX  D  -  PHASE  ONE:  SAMPLE  GEOMETRY 

TEST  RESULTS 


Sample  1:  n=4,  Triangular  Tab 


Converging  fracture  geometry 

Average  angle  of  convergence=6.8  degrees 

Effective  Fracture  Length=80mm  (Maximum) 


•  Converging  fracture  geometry 

•  Average  angle  of  convergence^  0.1 5  degrees 

•  Effective  Fracture  Length=41mm 
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. 


# 


# 


Sample  4:  n=6,  Trapezoidal  Tab 


•  Converging  fracture  geometry 

•  Average  angle  of  convergence=7.8  degrees 

•  Effective  Fracture  Length=80mm  (Maximum) 
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•  Converging  fracture  geometry 

•  Average  angle  of  convergence=8.25  degrees 

•  Effective  Fracture  Length = 80mm  (Maximum) 


m 
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Sample  Tensile  Test  Results 

0.71 1mm  Thickness  Sample 


0.406mm  Thickness  Sample 


APPENDIX  F  -  PHASE  TWO:  TEST  APPARATUS 
GEOMETRY  TEST  RESULTS 

Flush  Mounted  Geometry 


Sample  1:  Parallel  Cylinder ;  15mm  Radius 


Average  angle  of  convergerice=7.8  degrees 


•  Average  angle  of  convergence=5.65  degrees 
Sample  3:  Spherically  Tapered.  20mm  Maximum  Radius 


•  Average  angle  of  convergence-4.55  degrees 


# 
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Recess  Mounted  Geometry 


•  Average  angle  of  convergence=2.1  degrees 
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Sample  5:  Spherically  Tapered,  20mm  Maximum  Radius 


Average  angle  of  convergence=2.5  degrees 


APPENDIX  G  -  APPARATUS  DESIGN,  GEOMETRY  AND 

l 

SPECIFICATIONS 

Design  Details 
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Component  Specifications 

a  -  Top  Plate 

C-1018  Low  Carbon  Steel  -  lOinXIOinXO.Sin.  Conforms  to  ASTM  A108.  Heat- 
treating,  in  contact  with  carbon  (carburizing),  hardens  the  surface  of  this  low-carbon 
steel.  It  is  easy  to  cold  form,  bend,  braze,  and  weld.  Maximum  attainable  Rockwell 
hardness  is  B72.  Melting  point  is  2800  F.  Yield  strength  is  55,000psi.  Cold  finished. 
Width  and  length  tolerances  are  ±0.125in.  Thickness  tolerance  is  ±0.003in. 

b  -  Threaded  Adjustment  End 

Plain  Steel  Positive  Grip  Wire  Rope  End  Fittings.  Fitted  with  0.25in.,  28  thread, 
end  details. 

c  --  Wire  Tope 

Galvanized  Steel  Multi-Purpose  Rope — 7x19  class  strand 
core  commercial  grade.  Unlubricated  rope  offers  a  good  balance 
of  strength  and  flexibility  in  diameters  less  than  l/8in.  It  is  stronger 
but  less  flexible  than  six-strand  core  constructions. 

7  x  19  Class 
Strand  Core 

Galvanized  wire  rope  has  a  zinc  coating  that  provides  added 
corrosion  protection.  In  mild  environments,  it’s  an  economical  alternative  to  stainless 
steel.  The  strength  of  galvanized  rope  is  generally  less  than  that  of  plain  steel  and  stainless 
steel.  0.125in  diameter  and  20001b.  Breaking  strength.  Meets  specifications: 

•  Fed.  Spec.  RR-W-410 

•  Breaking  Strength  of  Mil-DTL-83420 

This  wire  rope  displays  the  following  linear  load-extension  characteristics: 

Load  =  Percent  Strain  *  Modulus 
where: 


Modulus— 1190.25  pounds  per  percent  extension 


d  -  Pillow  Block  Assembly 

C-1018  Low  Carbon  Steel  Precision  Ground  Stock  Conforms  to  ASTM  A108. 
Heat-treating,  in  contact  with  carbon  (carburizing),  hardens  the  surface  of  this  low- 
carbon  steel.  It  is  easy  to  cold  form,  bend,  braze,  and  weld.  Maximum  attainable 
Rockwell  hardness  is  B61-B62.  Melting  point  is  2800  F.  Yield  strength  is  55,000psi. 
Cold  finished.  Width  and  length  tolerances  are  +0.005in.  Thickness  tolerance  is 
±0.001in. 

e  -  Sample  Fastener  Bar 

f-  Double-Koto  Ball  Bearings 

Double-Row,  Double  Shielded 
Steel  Ball  Bearing  -  ABEC-1.  Double¬ 
row  ball  bearings  handle  high  radial  loads. 

The  balls  are  held  in  place  at  25°  angles 
between  the  inner  and  outer  sleeves.  They’re  ideal  for  pumps,  gear  motors,  and  large 
electric  motors.  Temperature  range  is  -40°  to  +250°  F.  Double-shielded  bearings  have 
steel  shields  that  help  keep  out  dirt  and  preserve  lubricants.  Shaft  diameter  20mm. 
Outside  Diameter  47mm.  Width  20.6mm.  Maximum  dynamic  radial  load  44501b. 
Maximum  RPM  10000. 

g  -  Tapered  Cylindrical  Koller 

12L14  Carbon  Steel  Rod.  Conforms  to  ASTM  A108.  Low-carbon  steel  that  has 
excellent  machining  characteristics  and  good  ductility  making  it  easy  to  bend,  crimp,  and 
rivet.  It  is  very  difficult  to  weld  and  cannot  be  case  hardened.  Maximum  attainable 
Rockwell  hardness  is  B75-B90.  Melting  point  is  2800°  F.  Yield  strength  is  60,000- 
80,000psi.  Cold  drawn. 

h  -  Base  Plate 

C-1018  Low  Carbon  Steel  -  12inX12inX0.5in.  Conforms  to  ASTM  A108.  Heat- 
treating,  in  contact  with  carbon  (carburizing),  hardens  the  surface  of  this  low-carbon 


steel.  It  is  easy  to  cold  form,  bend,  braze,  and  weld.  Maximum  attainable  Rockwell 
hardness  is  B72.  Melting  point  is  2800  F.  Yield  strength  is  55,000psi.  Cold  finished. 
Width  and  length  tolerances  are  ±0.125in.  Thickness  tolerance  is  ±0.003in. 


